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EXPERIMENTAL AND ANALYTICAL STUDY OF A CONICALLY DIFFUSED 
FLOW WITH A NEARLY SEPARATED BOUNDARY LAYER 
by Donald R. Boldman and Harvey E. Neumann 
Lewis Research Center 

SUMMARY 

Turbulence measurements were obtained in the nearly separated flow in a 13° total 
angle of divergence conical diffuser coupled to a constant area tailpipe. Air at 207 new- 
tons per square centimeter and 308 K provided an inlet velocity of about 51 meters per 
second at an inlet unit Reynolds number of 63.7x10° per meter. Very high longitudinal 
turbulence intensities accompanied the diffusion process with peak values approaching 
40 percent when normalized by the local centerline velocity. 

Comparison of the turbulence transport in the nearly separated boundary layer with 
values for flows with zero acceleration revealed that the ratio of turbulent shear to the 
turbulent kinetic energy was only about one-third of the value for boundary layers in 
zero acceleration. Also the peak values in the turbulent shear production and Reynolds 
stress generation terms were displaced further away from the wall in the case of the 
nearly separated boundary layer. 

Predictions of the pressxire recovery coefficient using a mixing length concept were 
good in the early stages of diffusion. Satisfactory predictions of the pressure recovery 
after predicted separation were obtained with an empirical constant effectiveness cri- 
terion. The best agreement was within 12 percent of the experimental pressure recov- 
ery. 


INTRODUCTION 

A knowledge of the turbulent transport associated with the subsonic diffusion of air 
is important in the design of components for high performance aircraft. With subsonic 
diffusion, a rise in static pressvire is obtained at the expense of a reduction in momen- 
tum. Eventually the flow in the decelerated boimdary layer near the wall has insufficient 
momentum to accommodate the increase in pressiire and the flow will tend to separate. 



During separation, viscous losses are high and consequently, the diffusion process be- 
comes very inefficient. In most practical diffusers, it is desirable to avoid separation; 
however, it is often necessary to operate the diffuser in a state of partial or incipient 
separation. The present study was performed in an effort to learn more about the tur- 
bulent structure of a diffusing flow which is in this state of incipient separation. 

In the present study a diffuser calculation method (ref. 1) based on the mixing length 
concept was applied to the flow in a conical diffuser having a nearly separated boundary 
layer. The predicted mean velocity profiles and pressure recovery are compared to the 
experimental results using a standard two- layer turbulence model. Results are also 
presented for a modified outer region turbulence model which was suggested on the 
basis of turbulence measurements in the diffuser. The constant effectiveness criterion 
Ox reference 2 is used to extend the pressure recovery prediction after separation. 

All of the present tests were conducted using a 13° total angle of divergence diffuser 
cbupled to a constant area tailpipe. The working fluid consisted of air at a nominal 
pressure of 207 newtons per square centimeter and a nominal total temperature of 308 K. 
A critical flow nozzle at the exit of the tailpipe was used to limit the inlet velocity to 
about 51 meters per second (inlet unit Reynolds number of 63. 7x10 /m) in order to stay 
within the incompressible flow regime. Hot-wire anemometers were used to measure 
the mean velocity profiles, longitudinal and lateral components of turbulence intensity, 
and the Re 3 molds stress. 


SYMBOLS 


A area 

a constant in hot-wire calibration equation 

,aj^ constant in hot-wire calibration equation 

b constant in hot-wire calibration equation 

bj^ constant in hot-wire calibration equation 

Cj skin friction coefficient 

Cp pressure recovery coefficient 

D diameter 

E bridge voltage 

e fluctuating component of voltage 

K constant in equation for the mixing length (eq. (7)) 
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constant in outer region turbulence model (eq. (9)) 

K 2 term in equation for mixing length (eqs. (7) and (8)) 

I mixing length 

M Mach number 

n mass fliix exponent in hot-wire calibration equation 
P pressime 

~2 ~2 ~2 2 

q turbulence kinetic energy function, q = u + v + w 

R radius 

Re Reynolds number 

T temperature 

U mean longitudinal velocity 

u longitudinal fluctuating component of velocity 

u’ root -mean -square (rms) value of u 

V effective velocity sensed by hot wire 

V lateral fluctuating component of velocity 

v' root-mean-square (rms) value of v 

w fluctuating component of velocity normal to x-y plane 

w' root-mean- square (rms) value of fluctuating component of velocity normal to u 

and v components 

X defined by eq. (BIO) 

X axial direction 

V defined by eq. (Bll) 

y distance from wall along normal 

a angle between wire number I of X-wire probe and mean flow direction 

/3 angle between wire number n of X-wire probe and mean flow direction 
y ratio of specific heats 

6 boundary layer thickness 

6* boundary layer displacement thickness 

e eddy viscosity 
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T 7 pressure recovery effectiveness, C /(C -j) 

P ' \ Pj IQ/j^ 

ju dynamic viscosity 

u kinematic viscosity 

p density 

a local half angle of divergence of diffuser 
9 boundary layer momentum thickness 

T shear stress 

(p angle between wires of X-wire probe and mean flow direction 
Subscripts; 

c compressible or centerline value 

e ambient conditions 

i incompressible 

id ideal 

in inner region of boundary layer 

I laminar 

o diffuser entrance 

out outer region of boundary layer 

ref reference value 

sep condition at separation 

t turbulent 

w condition associated with hot wire or wall 

9 based on momentum thickness 

0 stagnation condition 

Superscript: 

— average value 


APPARATUS 

The geometry, surface finish, and operating conditions for the present study were 
selected largely on the basis of the empirical design criteria of reference 3. In an 
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effort to obtain a nearly separated diffusing flow, a conical diffuser having a 13° total 
angle of divergence was chosen for this experiment. The diffuser was coupled to con- 
stant diameter pipe sections both upstream and downstream of the diffuser as shown in 
figure 1. The wall radii of curvature at the diffuser entrance and exit were equal to the 
entrance diameter D^, where = 7. 62 centimeters. The inlet, diffuser, and tailpipe 
were fabricated from AISI 304 stainless steel and were machined to a relatively smooth 
internal surface finish of 81.2 microcentimeters rms. 

Air entered the diffuser through a plenum, bellmouth, and 7.62-centimeter- 
diameter by 27. 3-centimeter-long pipe inlet as shown in figure 1. A3. 81-centimeter- 
diameter critical flow nozzle was located downstream of the tailpipe. Coordinates for 
the diffuser and tailpipe sections are given in table I. 



6l5° 


L27 



Hot-wire probe (typical, 5 places) 


Diffuser exit (wall detail) 


Diffuser entrance (wall detail) 



CO-11327-33 


Figure 1. - Apparatus for study of diffusing (low. Inlet diameter, 0^= 7.62 centimeters. (All dimensions in cm). 
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TABLE I. - DIFFUSER AND TAILPIPE 


COORDINATES 

(inlet diameter, D^ = 7.62 cm; inlet area, 
A^ = 45. 6 


station 

Axial 

distance, 

X, 

cm 

Diameter, 

D, 

cm 

Axial 

distance, 

x/Do 

Area 

ratio, 

A/Ao 

^0 

0 

7. 

620 

0 

1. 

000 

1 

4.209 

8. 

458 

.553 

1. 

232 

2 

6.749 

9. 

032 

.886 

1. 

405 

3 

9.296 

9. 

614 

1.220 

1. 

592 

4 

11.829 

10. 

191 

1.553 

1. 

789 

5 

14.369 

10. 

767 

1. 886 

1. 

997 

6 

16.911 

11. 

346 

2.220 

2. 

217 

7 

19. 451 

11. 

923 

2.553 

2. 

448 

8 

21.989 

12. 

497 

2. 886 

2. 

690 

^9 

24. 524 

13. 

071 

3.219 

2. 

942 

10 

27. 069 

13. 

653 

3.553 

3. 

210 

11 

29. 609 

14. 

229 

3.886 

3. 

487 

12 

32. 151 

14. 

803 

4.220 

3. 

774 

13 

34.686 

15. 

380 

4.552 

4. 

074 

14 

37.229 

15. 

956 

4. 886 

4. 

385 

15 

39.769 

16. 

528 

5.219 

4. 

705 

16 

42.314 

17. 

104 

5.553 

5. 

039 

^17 

44. 869 

17. 

686 

5.889 

5. 

387 

(b) 

54. 864 

19. 

533 

7.200 

6. 

571 

18 

62.484 



8.200 



^19 

77. 089 



10. 117 



20 

91.694 



12.034 



21 

106.299 



13.950 



^22 

120. 904 



15.867 




^Denotes boundary layer survey station. 
'’Tange ncy. 


INSTRUMENTATION 
Wall Static Pressures 

Pressure taps having a diameter of 0. 0794 centimeter were located in the diffuser 
and tailpipe starting at the diffuser entrance (x/D^ = 0 in fig. 1). A pressure tap was 
located at each of the twenty-three stations listed in table I. 
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The static pressures were measured by means of manometers containing acetylene 
tetrabromide which has a nominal room temperature density of 2. 934 grams per cubic 
centimeter. The stagnation pressure, measured at the plenum centerline, provided the 
reference pressure for the manometer system. 


Boundary Layer Surveys 

Hot-wire measurements of mean velocity, longitudinal and lateral components of 
turbulence intensity, and Reynolds stress were obtained at stations 9 and 17 in the dif- 
fuser and at stations 19 and 22 in the tailpipe (refer to table I). A single wire probe was 
used for the measurements of mean velocity and longitudinal components of turbulence 
intensity, whereas an X-wire probe was used to measure the lateral components of tur- 
bulence intensity and the Reynolds stress. Since the boundary layer at the diffuser en- 
trance was too thin to accommodate the X-wire probe, only the mean velocity and longi- 
tudinal turbulence intensity distributions were measured at this location. 

All of the hot-wire measurements were obtained with constant temperature hot-wire 
anemometers using 0. 0005 -centimeter- diameter tungsten hot wires. The hot-wire 
probes were moved point -by- point from the centerline to the wall by motorized actuators. 
The single wire probes were driven until contact with the wall was established. The X- 
wire probe traverse was terminated prior to contact with the wall by means of a sleeve 
located on the probe shaft. The minimum distance from the wire crossing point to the 
wall was 0.381 centimeter. Further details concerning the actuators and single wire 
probe design can be obtained from reference 4. 

A pitot pressure probe having a flattened tip 0. 005 centimeter high by 0. 076 centi- 
meter wide was also used to measure the velocity profile at x/D^ =3.22 in the diffuser. 
Details concerning the pressure probe design can be obtained from reference 5. 


EXPERIMENTAL FLOW CONDITIONS 

All data were obtained with the same nominal flow conditions. The stagnation pres- 
sure Pg was maintained at 207±1 newtons per square centimeter absolute. The stag- 
nation temperature Tq for the tests was 308±3 K. These stagnation conditions provided 

a nominal inlet velocity of 51.3 meters per second and an inlet Resmolds number 

fi ^ 

based on D of 4. 85x10 . 

0 
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EXPERIMENTAL RESULTS 


Pressure Recovery 


The experimental pressure recovery in the diffuser and tailpipe is usually presented 
in terms of the pressure recovery coefficient Cp where 



( 1 ) 


The velocity is the mass-averaged value at the diffuser inlet. The experimental 
values of Cp are often compared to the maximum attainable values for a one- 
dimensional flow of a perfect gas in a diffuser of the same geometry. These ideal values 
of C can be expressed as 



( 2 ) 


for incompressible flow and as 






-y/(y-l) 



y- 1 
2 



-y/(y-l) 





-y/(y-l) 


( 3 ) 


for compressible flow. The experimental and ideal distributions of pressure recovery 
coefficient for the diffuser and tailpipe are presented in figure 2 and are tabulated in 
table II. The distributions of ideal pressure recovery, (Cp and (Cp , in fig- 


ure 2 indicate that the flow in this investigation can, for all practical purposes, be con- 
sidered as incompressible. A maximum ideal pressure recovery coefficient of 0.977 is 
obtained at the diffuser exit and remains constant in the tailpipe (refer to eq. (2)). The 
experimental pressure recovery coefficients are appreciably lower than the ideal values 
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with Cp = 0. 72 at the diffuser exit (x/D^ = 7. 2). The tailpipe provided only a modest 
increase in Cp of from 0. 72 at the diffuser exit to about 0. 75 near the tailpipe exit. 

It is recognized that in the entrance region of the diffuser there is a modest favor- 
able pressure gradient which can influence the value of (P^/Pq) used as the reference 

pressure in the pressure recovery coefficient (eq. (1)). An indication of the magnitude 
of this favorable pressure gradient in the present diffuser will be given in a later sec- 
tion entitled PREDICTED RESULTS. 

The solid symbols in figure 2 represent the values of Cp at the boundary layer 
survey stations. It should be noted that a boundary layer survey was obtained at the 
diffuser entrance in the plane of the static pressure tap which was used for the reference 
pressure in computii^ Cp. The mean velocity profile at this station provided one of the 
input matching conditions for the diffuser flow analysis of reference 1. 


Experiment 

Ideal, compressible; y ■= 1.4 


Ideal, incompressible 

Solid symbols denote boundary layer survey stations 



Figure Z - Pressure distribution in diffuser and tailpipe. Inlet Reynolds 
number, Re^ = 63. 7x10^ per meter; inlet diameter. Dp = 7. 62 centimeters. 
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TABLE n. - EXPERIMENTAL AND IDEAL PRESSURE RECOVERY COEFFICIENTS 


[inlet diameter, D^ = 7.62 cm; stagnation pressure, Pq = 207 N/cm^O 


Station 

Axial 

distance, 

x/Dq 

Pressure 

ratio, 

Pw/Po 

Mach 

number, 

M 

Pressure recovery coefficient 

Experimental 

S 

Ideal incom- 
pressible 

S, id i 

Ideal incom- 
pressible 

S.idc 

^0 

0 

0.98528 

0.146 

0 

0 


0 


1 

.553 

.98985 

. 121 

.304 

.341 

.347 

2 

. 886 

.99152 

.110 

. 416 

.493 

.500 

3 

1.220 

.99253 

. 104 

.484 

.605 

.613 

4 

1. 553 

.99323 

.0986 

.531 

.687 

.695 

5 

1.886 

.99374 

.0948 

.566 

.749 

.756 

6 

2.220 

. 99406 

.0923 

.589 

.797 

.804 

7 

2.553 

.99440 

.0896 

.611 

.833 

.840 

8 

2.886 

.99464 

.0877 

.627 

.862 

.868 

^9 

3.219 

. 99487 

.0858 

.643 

.885 

.891 

10 

3.553 

. 99503 

.0844 

.654 

.903 

.909 

11 

3. 886 

.99529 

.0833 

.663 

.918 

.924 

12 

4.220 

. 99529 

.0822 

.672 

.930 

.936 

13 

4. 552 

. 99542 

.0810 

.680 

.940 

.946 

14 

4.886 

. 99554 

.0799 

.689 

.948 

.954 

15 

5.219 

.99561 

.0793 

.693 

.955 

.961 

16 

5.553 

.99568 

.0787 

.698 

.961 

.966 

^17 

5.889 

. 99577 

.0778 

.704 

.966 

.971 

18 

8.200 

.99617 

.0741 

.731 

.977 

.983 

^19 

10. 117 

. 99637 

.0721 

. 745 





20 

12.034 

. 99646 

.0712 

.751 





21 

13.950 

.99651 

.0708 

.754 





^22 

15.867 

.99651 

.0708 

.754 






^Denotes boundary layer survey station. 



Boundary Layer Surveys 


In this section, the results of hot-wire surveys of the' mean velocity U, longitudinal 
component of turbulence intensity u', lateral component of turbulence intensity v', and 
the Reynolds stress uv will be presented. The turbulence measurements were ob- 
tained with single and X-Avire probes. A linearized circuit was used in conjunction with 
the single wire measurements to obtain the mean velocity profiles. Details concerning 
the hot-wire data reduction are presented in appendix A. As in the diffuser experiment 
of reference 6, the accuracy of the turbulence measurements was difficult to establish. 
In reference 6, the estimated errors which were based on the results of reference 7 
were 5 percent in u' and 10 to 15 percent in uv during the early stages of diffusion. 

In the latter stages of diffusion the estimated error increased to about 20 percent in u' 
and 20 to 40 percent in uv. The latter stage of the diffusion process in the conical dif- 
fuser of this study probably occurs downstream of station 9 (x/D^ = 3.22). The bound- 
ary layer data will be presented graphically but will not be tabulated because of the dif- 
ficulty in determining the measurement error and the possibility of establishing a false 
sense of the measurement accuracy. 

Another source of error in the turbulence measurements arises because of mis- 
alinement of the X-wire probe which results in unequal angles between the wires and the 
mean flow velocity vector. This error is present at station 9 (x/D^ =3.22) and at sta- 
tion 17 (x/Dq = 5.89) in the diffuser because the probe is installed in a position normal 
to the wall. The effect of this misalinement error is to reduce the measured values of 
v' and uv with the maximum reduction occurring at the centerline. The effects of this 
misalinement error will be ignored in the presentation of the turbulence data; however, 
the magnitude of the required correction at the centerline of stations 9 and 17 has been 
estimated in appendix B. 

Before examining the profiles of mean velocity and turbulence intensity, consider 
the flow development in the diffuser and tailpipe in a qualitative manner. The flow in 
the present system is depicted schematically in figure 3. Actually, this figure is based 
on the profile data yet to be presented. The large inviscid core in the entrance region 
(x/D^ = 0) rapidly becomes displaced by a turbulent core. The complete disappearance 
of the inviscid core takes place in the region bounded by 3. 22 < x/D^ < 5. 89 - probably 
just slightly downstream of x/D^ = 3. 22. In analyzing the results, emphasis will be 
placed on the data at x/D^ = 3. 22 since this station is the first one to reflect the 
effects of the strong adverse gradient in the entrance region. 

Velocity profiles . - The boundary layer velocity profiles at the three stations in the 
diffuser and at the two stations in the tailpipe are presented in figure 4. The boundary 
layer velocity U is nondimensionalized by the local centerline value U^(x) which is 
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Axial 



Choked nozzle 



Figure 4. - Mean velocity profiles in diffuser and tailpipe. 
Inlet Reynolds number, Re^^ 63. 7x10^ per meter; inlet 
diameter, Dg = 7.62 centimeters. 


also the maximum velocity at a given station. The value of Up(x) decreases from 
51.3 meters per second at the entrance to 7.7 meters per second at the downstream 
station in the tailpipe. The local radius R(x), which is used to nondimensionalize the 
distance variable y, increases from 3. 81 centimeters at the entrance to a maximum of 
9.77 centimeters in the tailpipe. 
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Referring to figure 4, it can be noted that the boundary layer at the- diffuser en- 
trance is relatively thin, 6 r; 0. 2 R(x); however, as the flow diffuses the boundary layer 
rapidly fills the channel. This rapid growth can be noted by comparing the profiles at 
x/Dq = 0 and 3.22. The profiles at x/D^ = 3. 22 and 5. 89 suggest that the boundary 
layer is in a nearly separated state. As in reference 8, the velocity near the wall ap- 
proaches a value of about 20 percent of the local maximum U(.(x). Furthermore the 
boundary layer profile near the wall is similar to the classical incipiently separated 
shape, that is, 9U/3y - 0. Further development of the viscous flow in the tailpipe re- 
sults in a readjustment or relaxation to a nearly uniform profile across the duct 
(x/Djj = 15.87). 

A pitot pressure probe was also used to survey the velocities at x/D^ = 3. 22. Al- 
though the shape of the profile is different very close to the wall (y/R(x) < 0. 05), the 
pitot probe results are generally in agreement with the hot-wire measurements and indi- 
cate the presence of incipient separation rather than actual separation. 

Longitudinal component of turbulence intensity . - Radial distributions of the longi- 
tudinal component of turbulence intensity u' are represented in two different nondimen- 
sional forms, u'/U^(x) and u'/U^, and are presented in figures 5 and 6, respectively. 



Figure 5. - Distributions of longitudinal component of turbulence intensity 
(normalized by local centerline vetocity) in diffuser and tailpipe. Inlet 
Reynolds number, Re^^ 63. 7x10^ per meter; inlet diameter, Oq- 7.62 
centimeters. 
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Axial Radius, 



Diffuser 

centerline 


Figure 6. - Distributions of longitudinal component of turbulence intensity 
(normalized by entrant velocity) in diffuser and tailpipe. Inlet Reynolds 
number, Reg =63. 7xlOP per meter; inlet diameter, Dg = 7. 52 centimeters; 
inlet velocity, Ug = 51.3 meters per second. 


The nondimensionalization by provides a representation of the change in the abso- 
lute value of u’ since U^, the diffuser entrance velocity, is a constant whereas the 
nondimensionalization by U^(x) gives an indication of the radial change in u' relative 
to the local centerline velocity. The same nondimensionalization will also be used in 
presenting the results of the lateral components of turbulence intensity v' and the 
Re37nolds stress uv. 

The radial distributions of the longitudinal component turbulence intensity u'/U^(x) 
at the five stations are shown in figure 5. The profile at the diffuser entrance indicates 
a free stream intensity of about 1. 0 percent and a maximum near the wall of about 
8. 0 percent. As the flow diffuses, the intensity increases appreciably with peak values 
of u'/U^(x) approaching 40 percent (x/D^ = 10. 12). The maximum in the intensity dis- 
tribution moves toward the centerline as the flow diffuses. Near the exit of the tailpipe 
where the flow appears to be reestablishing some form of equilibrium, the maximum 
value of u'/U^(x) occurs near the wall. This trend of increasing turbulence intensity 
with decreasing axial velocity has also been observed by Spangenberg, Rowland, and 
Mease (ref. 9) in a two-dimensional diffusing flow near separation. In reference 9, 
values of u'/U in excess of 0. 5 were obtained in a region very near the wall. Similar 
levels of u'/U occur in the wall region of the present diffuser at x/D^=3.22. How- 
ever, during the latter stages of diffusion, that is, at x/D^ = 5.89 and 10. 12, these 
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high levels of intensity (u'/U > 0. 4) extend from the wall region to approximately the 
duct centerline. This can be noted by dividing the intensity u’/U^(x) in figure 5 by the 
corresponding velocity ratio U/U^(x) in figure 4. 

The distribution of u'/U^ in figure 6 indicate that the absolute value of u' also 
increases as the flow passes from the entrance to the diffuser exit (x/D^ = 5. 88). The 
measurements indicate that the maximum value of u’ is about 16 percent of the en- 
trance velocity or about a factor of two higher than the observed maximum in u' at the 
diffuser entrance. This change in u'/U is much less than the change in u'/U (x) 

(fig. 5), thus indicating that a large part of the apparent increase in the intensity 
u'/Up(x) is the result of nondimensionalization by the local centerline velocity. A simi- 
lar effect involving the method of nondimensionalizing u' has beep observed in the 
boundary layer of accelerating flows (e.g. , ref. 10). 

Lateral component of turbulence intensity . - The radial distributions of the lateral 
component of turbulence intensity v' are presented in the two nondimensional forms in 
figures 7 and 8. Comparison of these results with the corresponding values of u' in 



Distance from wall. y/R(x) 


Figure 7. - Distributions of transverse component of turbulence intensity 
(normalized by local centerline velocity) in diffuser and tailpipe. Inlet 
Reynolds number, ReQ = 63. 7x10^ per meter; inlet diameter, Dq = 7.62 
centimeters. 
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Distance 

centerline 


Figure 8. - Distributions of transverse component of turbulence intensity 
(normalized by entrance velocity) in diffuser and tailpipe. Inlet Reynolds 
number, Re^ =63.7x10^ per meter; inlet diameter, Dq • 7. 62 centimeters; 
inlet velocity, Ug ■ 51. 3 meters per second. 


figures 5 and 6 indicates that v' is much less than u'. In the latter stages of diffusion, 
the maximum in v' moves towards the centerline. The maximum values of u' and v' 
at a given station occurred at about the same radial position. 

Reynolds stress distributions . - The distributions of the Reynolds or shear stress 

functions uv/Uj.(x)^ and are presented in figures 9 and 10, respectively. A 

maximum value of uv/u^(x)^ of about 0.0085 occurs at x/D^ = 10. 12 (fig. 9) whereas 
the maximum in the absolute value of uv occurs at x/D^ =3.22 as shown in figure 10. 
These values of uv must be adjusted, however, for the errors mentioned previously 
and discussed in appendix B. In the next section the distribution of uv at x/D^ =3.22 
will be used in conjunction with the corresponding mean velocity profile to estimate the 
constant in the outer region of the Clauser (ref. 11) eddy viscosity model used in the 
analysis of reference 1. 
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centimeters diameter, Dq = 7. 62 centimeters; inlet velocity. Uq - 51.3 

meters per second. 



Estimation of Eddy Viscosity at x/D^ = 3.22 

It is well known that calculations of the turbulent boundary layer require a certain 
amount of empirical information concerning the turbulence modeling. Upon specifpng 
the turbulence transport, the boundary layer partial differential equations of motion can 
be integrated directly by numerical methods. The turbulence model of interest in this 
study is based on the mixing length concept in which the turbulence or eddy viscosity is 
related to the mean flow properties. In the boundary layer analysis of reference 1, 
which has been applied to the present data, the eddy viscosity in the boundary layer is 
described by a two- layer model. This model will be reviewed briefly because of its 
importance relative to the predicted boundary layer in the present diffuser. 

The shear stress term appearing in the equations of motion can be expressed as the 
sum of a viscous and Reynolds stress contribution; namely, 

T = T, + Ti. = M — - PUV (4) 

^ ^ ay 


where 


-puv = pe — (5) 

ay 

The eddy viscosity e is related to the Reynolds stress uv and the mean velocity gradi- 
ent au/ay as shown in ec[uation (5). 

The empirical model for the eddy viscosity appearing in the analysis of reference 1 
is based on a two-layer concept. In the two- layer model, the eddy viscosity in the near 
wall region is expressed as 


€ 



9y 


where the mixing length I is given by the VanDriest (ref. 12) expression: 


I 




( 6 ) 


(7) 
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The value of K 2 in equation (7) is given by 


/ 

K2=26(^) (8) 

W/ 


The Clauser (ref. 11) outer region model used in the analysis of reference 1 can be ex- 
pressed as 


e 


^out 


1 c 


( 9 ) 


The constants K and in equations (7) and (9) are 


K = 0. 41 ^ 

= 0. oiej 


( 10 ) 


The mixing length model for the inner portion of the boundary layer depends on the 
wall shear stress which was not measured in the present study. Therefore a direct 
comparison of e based on equations (6) and (7) with the values of e determined from 
the measured velocity and Reynolds stress distributions cannot be made; however, this 
is not the case for the outer region of the layer. As shown in equation (9), the eddy 
viscosity is based on the mainstream velocity and the displacement thickness 6*, 
where 


6 * 



U(x, y) 
U^(x)_ 



—S— cos a 
R(x) J 


dy 


( 11 ) 


The eddy viscosity function e/6* U^(x) for the outer region of the boundary layer 
was determined from the boundary layer measurements at x/D^ = 3.22. This station 
was selected because it represents the only station downstream of the diffuser entrance 
in which the boundary layer has not completely filled the channel (refer to fig. 3). The 
results of e/6* Uj.(x) are plotted as a function of y/R(x) in figure 11. The standard 
constant Kj (Kj = 0. 016) in the analysis of reference 1 is shown as a dashed line. 

Although there are inaccuracies in the measured values of uv, an eddy viscosity 
distribution can be presented. Clearly, the experimental distribution of Kj in figure 11 
differs from the constant value of 0. 016 which was used in the Clauser outer region 
model (ref. 11). The Clauser model does not appear to apply to the present results be- 
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Figure 11. - Eddy viscosity distribution based on measured mean velocity 
and Reynolds stress distributions at axial distance x/D^ = 3. 22. Bound- 
ary layer displacement thickness, 5* = 1.45 centimeters; centerline ve- 
locity, U(.(x) = 30. 2 meters per second; radius, R(x) -6.54 centimeters. 


cause it underestimates the outer region turbulence transport. Alternatively, in the 
absence of a more definitive model, the results of figure 11 suggest that the outer region 
constant Kj might be increased to account for the increased turbulence. The effect 
of Kj on the predicted pressure recovery coefficients and boundary layer velocity 
profile at x/D^ = 3. 22 will be presented in the following section of the report. 


PREDICTED RESULTS 

It is recognized that the unsteady effects associated with incipient separation are 
not taken into account by conventional equilibrium turbulence models. Although the 
neglect of the unsteady phenomena is valid for a large number of nonzero pressure gra- 
dient flows, this may not be the case for flows at incipient separation. Unfortunately 
these unsteady effects are not well understood and therefore in order to estimate the 
boundary layer in such flows it is necessary to resort to existing turbulence models. 

With this in mind, the sophisticated finite-difference diffuser calculational method of 
reference 1 has been applied to the present flow in an effort to predict the pressure re- 
covery and mean velocity profiles. 

The method of reference 1 solves for the entire flow across the duct at each stream- 
wise station. The problems of matching solutions for the inviscid flow and boundary 
layer are thereby eliminated and no difficulties exist when the boundary layer merges 
at the diffuser centerline. This analysis assumes that the streamlines for the actual 
flow through the duct will not be appreciably different from the potential flow stream- 
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lines, and thus it is possible to make boundary layer-type approximations in a coordi- 
nate system based upon the stream function and velocity potential of the potential flow. 
The viscous effects are therefore treated as a perturbation upon the inviscid flow field 
and are governed by parabolic partial differential equations. 

Experience with the program indicated that the wall curvature at the diffuser en- 
trance caused numerical difficulties, when the calculation was initialized at station 0 
(table I). In applying the analysis to the present configuration, it was convenient to 
initialize the calculations upstream of the entrance, that is, in the pipe section. By 
starting the calculation in this way the numerical problems in the entrance were avoided 
and the calculation provided an estimate of the real flow conditions at the diffuser en- 
trance. 


Pressure Recovery Coefficient 

Entrance region . - The results of the analysis of reference 1 which are based on a 
reference pressure in the pipe (x/D^ = - 2. 59), indicate a negative value of up to the 
diffuser entrance as shown in figure 12. The difference between the experimental and 
theoretical values of Cp at the diffuser entrance results from the approximations used 
in initializing the calcuMion. The error in Cp resulting from the uncertainty in the 
inlet pressure should not exceed 0. 04. 



Figure 12. - Predicted pressure distribution in diffuser en- 
trance region. Iniet Reynolds number, ReQ = 63.7xlO^ 
per meter; inlet diameter, Dq = 7. 62 centimeters. 
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Figure 13. - Comparison of predided and experimental pressure recovery coefficients 
in diffuser for various values of Kj in outer region eddy viscosity model. Inlet 
Reynolds number, Re,, = 63. 7x# per meter; inlet diameter. Dp = 7. 62 centimeters. 


Diffuser . - Predictions of the pressure recovery coefficient in the diffuser are com- 
pared to the experimental results in figure 13. The effect of increasing the constant 
in the outer region eddy viscosity model (eq. (9)) was to extend the predicted separation 
point further downstream and to increase the predicted values of C . Predicted sepa- 
ration occurred upstream of the survey station at x/D^ = 3.22 when Kj was equal to 
0. 010 and 0. 016. As noted in figure 13, the calculation of Cp in the upstream portion 
of the diffuser was insensitive to the value of Kj. In the early stages of diffusion the 
best agreement with experimental values of Cp was obtained with the low value of 
(Kj = 0. 010); however, separation was predicted very early in the expansion process. 
Use of the large value of Kj (Kj = 0. 030) resulted in a pronounced extension of the un- 
separated flow region; however, the predicted values of C in this extended region 
were appreciably higher than experiment. At the predicted separation point, C ex- 
ceeds the experimental value by about 18 percent. Although the larger value of 
appears more consistent with the data, the modeling of the outer eddy viscosity region 
exclusive of the inner region is not sufficient to accurately calculate through the nearly 
separated flow regime. It will be shown that estimates of C in this flow regime can 
be made with reasonable accuracy by means of the empirical constant effectiveness cri- 
terion suggested in reference 2 . The results based on this simple method are influenced 
by the location of predicted separation. In diffusers of the present type, where incipient 
separation rather than actual separation occurs, the calculated separation point must be 
treated hypothetically since it merely delineates a region in which the turbulence model- 
ing in the analysis breaks down. Further evidence of these deficiencies in turbulence 
modeling can be noted by examining the velocity profile development in the diffuser. 
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Velocity Profiles 


The predicted velocity profiles between x/D^ = 0 and 3.22 (station 9) are shown in 
figure 14. As indicated previously, in order to extend the calculation to x/D^ = 3. 22 
it was necessary to use large values of Kj (Kj = 0. 030). The calculated velocity profile 
at x/Dq =3.22 tends to attain the same general shape as the profile at incipient separa- 
tion; however, the local velocities are appreciably higher than the experimental values. 
These higher local velocities imply a thinner boundary layer and higher C , consistent 
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Figure 14. - Predicted boundry layer development in diffuser for different 
values of Kj in outer region edcfy viscosity model. Note that separa- 
tion was predicted upstream of boundary layer survey at x/Dj, = 3. 2 
when K j = 0. 010 and 0. 015. 
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with the results of figure 13. The principal effect of increasing Kj is to increase the 
local velocity near the wall and, consequently, delay predicted separation. 


Calculation of Cp After Predicted Separation (Ref. 2) 

In reference 2 it was shown that after the predicted separation point in a conical 
diffuser the effectiveness tj remains nearly constant, where 



In applying this result to the present data it was necessary to determine tj 

Ow U 

and simply multiply this constant by the ideal pressure recovery coefficient (C 


= constant 
p,id).- 


The results from this method are given in figure 15. The best agreement was obtained 


when separation was predicted at x/D^ =1.8 corresponding to the calculation with 
Kj = 0. 010. In this case the agreement was within about 12 percent of the experimental 
values of Cp. The largest discrepancy occurred at x/D^ 5 which is near the dif- 
fuser exit. Predicted results for Kj = 0. 016 and 0. 030 are about 6 percent higher than 
the predictions based on Kj = 0. 010, or within about 18 percent of the experimental 
values of Cp. 



Axial distance from diffuser entrance, x/Dq 


Figure 15. - Predicted pressure recovery coefficients after separation based 
on the constant effectiveness method (ref. 2). Re^, =63. 7x#/m; = 

7.62 centimeters. Shaded symbols denote boundary layer survey stations. 
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COMPARISON OF TURBULENCE TRANSPORT IN ZERO PRESSURE 


GRADIENT AND DIFFUSED BOUNDARY LAYERS 


In the presentation of the results of the turbulent shear stress measurements em- 
phasis was placed on the outer wake region of the boundary layer. It is of interest to 
further compare the turbulence results for the nearly separated boundary layer at sta- 
tion 9 (x/Dq 3.22) with the turbulence characteristics for a zero pressure gradient 
boundary layer. In making this comparison, the distributions of turbulence transport 
for the nearly separated boundary layer will be compared to the results for a zero pres- 
sure gradient boundary layer of thickness 6 = R(x) using the turbulence measurements 
of Klebanoff (ref. 13). The results of reference 13 were used as input in the data re- 
duction program of the present study in order to provide a consistent method of calcu- 
lating the transport terms for the two boundary layers. The thickness of the boundary 
layer at x/D^ = 3. 22 is considered to be equal to the duct radius R(x). 

The transport terms to be considered are as follows: -uv^^ q^, -uv R(x)(p/t^)^^^ 

[au(x)/3y], R(x)(p/t^)^^^ [au(x)/3y], and -uv/ v^. These terms represent respec- 

tively, (1) the ratio of turbulent shear stress to the turbulent kinetic energy, (2) the 
production of turbulence, (3) the generation of Reynolds or shear stress, and (4) the 
ratio of the production of turbulence to the generation of Reynolds stress. Since the w' 
component of turbulence intensity and the wall shear stress were not measured, it 
was necessary to make certain assumptions in order to obtain these quantities. The 
turbulent kinetic energy is 

q^ = u^ + v^ + w^ (13) 

2 

where w is obtained from the approximation that 


_ (u» + V*) 
2 


(14a) 


or 


„2 (u^ + 2uV + v^) 

w = - 

4 


(14b) 
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This assumption is quite good for the turbulent boundary layer in a zero pressure gradi- 
ent (see ref. 13) and in a mild adverse pressure gradient (ref. 14). The wall shear 
stress was calculated for convenience from the flat plate friction law in which 

Cj = 0.012 Reg°'^^ (15) 


and 


T... = — pUp(x) 
2 ^ 


w 


(16) 


The boundary layer momentum thickness 6 was given by 


e = 


/»R(x) 



_ 

/ U(x, y) 

1 _ 

U(x,y) 

1 - y rr 

i 

X 

U,(x)_ 

. R(x) 


dy 


(17) 


The assumption that the shear stress can be calculated from equations (15) and (16) is 
questionable since the boundary layer is in a nearly separated state. In the case of 
complete separation Cj would be zero; however, it is obvious from equations (15) and 
(17) that a nonzero value of Cj would be estimated. The calculated value of in the 
turbulence and Reynolds stress production terms is probably too high and therefore 
would tend to suppress the magnitude of these terms. 

The distribution of -uv/ q^ is presented in figure 16. Bradshaw (ref. 15), Hinze 

— l~2 

(ref. 16), and Townsend (ref. 17) indicate that the distribution of -uv/ q is nearly 

uniform over the large wake portion of the boundary layer and that -uv^ q^ attains a 
value of about 0. 14 to 0. 15. These values are consistent with the calculated results 
using Klebanoff's data (ref. 12) for the nonaccelerated boundary layer. The greatly re- 
duced values of -iT^ q^ for the diffuser boimdary layer imply that treating the viscous 
development in the outer portion of the boundary layer as a nonaccelerated boundary 
layer wake region (as in the analysis of ref. 1) may not be entirely adequate. A constant 
ratio of shear stress to kinetic energy also implies a similarity between the two distri- 
butions consistent with Von Karman's similarity hypothesis concerning the structure of 

turbulence (ref. 16). Although the level of the distribution of -u^ q^ for the diffuser 

flow is much lower than the zero acceleration case, the value of -u^q^ is approxi- 
mately constant over a large portion of the diffuser boundary layer (values ranging from 
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Figure 16. - Comparison of distributions of for 

nonaccelerated and diffused boundary layers. Inlet Reyn- 
olds number, Rep = 53.7x10^ per meter; inlet diameter, 

Dp = 7. 62 centimeters; radius, R(x) = 6. 54 centimeters. 


about 0. 04 to 0. 05). 

The distribution of the terms representing the production of turbulence and Reynolds 
stress are shown in figures 17 and 18, respectively. In the nonaccelerated boundary 
layer these production terms reach a peak near the wall whereas in the diffused bound- 
ary layer near separation the maximum values for the production terms occur much 
further from the wall (y/R(x) « 0. 4). It is also apparent from figures 17 and 18 that 
there is a strong viscous action in the diffused boundary layer which results in the pro- 
duction of turbulent shear over an appreciable extent of the channel. 

The ratio of the distributions of figures 17 and 18, -u^ v^, is plotted in figure 19. 
Note that any ambiguity in the production terms resulting primarily from the uncertainty 

in has been removed in forming the ratio of -u^ v^. In the nonaccelerated bound- 
ary layer the ratio of turbulence to Reynolds stress production is nearly constant over 
a large portion of the boundary layer; however, in the diffused boundary layer a distinct 

maximum occurs at y/R(x) « 0. 2. In general, the ratio of is lower in the 

diffused boimdary layer, especially in the outer portion of the layer. 

Examination of figures 16 and 19 reveals that there is a significant change in the 
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Figure 17. - Comparison of distributions of -uvRIx)^:^^ 

for nonaccelerated and diffused boundary layers. Inlet Reynolds 
number, Rej,=:63.7xlC^p€rmeter! inlet diameter, Dp'7.62 
centimeters; radius, R(x) “ 6. 54 centimeters. 



0 .2 .4 .6 .8 1.0 


Distance from wall, y/Rix) 

Figure 18. - Comparison of distributions of v^RIxl^:^^^ 

for nonaccelerated and diffused boundary layers. Inlet Reynolds 
number, Reg==63.7xl0^per meter; inlet diameter, Dq = 7.62 
centimeters; radius, R(x) = 6. 54 centimeters. 
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Figure 19. - Comparison of distributions of -uv /v^ for 
nonaccelerated and diffused boundary layers. Inlet Reyn- 
olds number, Reg =^63.7 x 106 per meter; inlet diameter, 

D() = 7. 62 centimeters; radius, R(x) = 6. 54 centimeters. 


turbulence structure from that associated with a zero pressure gradient flow. The low 
/~2 

values of -uv/ q for the adverse pressure gradient imply a relatively low turbulent 

' 12 

shear for a given turbulence intensity level. The turbulent shear term, -uv/ q is only 

about one-third of the value for a zero pressxire gradient flow (fig. 16). The transport 

of Reynolds stress by the mean flow depends on the generation as given by the product of 

~2 

the mean velocity gradient 3U/0y and the transverse turbulence component v (ref. 15). 
Likewise, the transport of turbulent kinetic energy by the mean flow depends on the 
production of turbulence as given by the product of the Reynolds stress and the mean 
velocity gradient. Therefore the low value of the ratio of turbulence production to shear 

production -u^ v^ contributes to a lower change in turbulent energy relative to the 
change in Resmolds stress. The ratio of these production terms is generally from about 
0.3 to 0.8 of the value for zero pressure gradient flow (fig. 19). Additional measure- 
ments would have to be made before conclusions can be drawn concerning the relative 
balance between the generation of turbulence and shear and the balance with dissipation, 
advection, and turbulent diffusion. 


SUMMARY OF RESULTS 

An experimental study was performed to determine some of the turbulence charac- 
teristics of the nearly separated flow in a 13° total angle of divergence conical diffuser 
coupled to a constant area tailpipe. The working fluid was air at a nominal pressure of 
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207 newtons per square centimeter and a nominal total temperatxire of 308 K. A critical 
flow nozzle at the exit of the tailpipe was used to limit the inlet velocity to about 51 me- 
ters per second (inlet unit Reynolds number of 63. 7x10®/“)* The experimental bound- 
ary layer velocity profiles and pressure recovery coefficients were compared to theo- 
retical results based on a diffuser boundary layer method which employs the mixing 
length concept. An empirical method called the constant effectiveness criterion was 
used to extend the prediction of pressure recovery downstream of the predicted separa- 
tion point. The results of this study can be summarized as follows: 

1. Hot-wire measurements of the velocity profiles at two stations in the diffuser 
indicated that the boundary layer was in a nearly separated state. This nearly separated 
condition persisted into the tailpipe recovery section. 

2 . The diffusion process was accompanied by extremely high values of longitudinal 
turbulence intensity with uVU^(x) approaching 40 percent in the latter stages of diffu- 
sion. The maximum value of u'/U^ at the diffuser entrance was about 8 percent and 
occurred very close to the wall; however, as the flow diffused, the maximum value of 

u'/U„(x) moved towards the centerline. 

c 

3. The ratio of turbulent shear stress to turbulent kinetic energy uv/q was only 
about 30 percent of the accepted values for a boundary layer with a zero or mild adverse 
pressure gradient. 

4. The maximum values for the tvirbulence production and Reynolds stress produc- 
tion for the nearly separated botmdary layer occurred at an appreciable distance from 
the wall relative to the position of the maximum for a nonaccelerated boundary layer. 

5. In the upstream portion of the diffuser good predictions of the pressure recovery 
coefficient were obtained with the diffuser boundary layer program; however, it is be- 
lieved that separation was predicted prematurely. The value of the constant in the outer 
region eddy viscosity model was increased from the standard value of 0. 016 to 0. 030 to 
be more consistent with the observed turbulence structure. The effect of using a higher 
value of the constant was to extend the calculation an appreciably greater distance down- 
stream but at the expense of pressure recovery accuracy. The resulting pressure re- 
covery was about 18 percent higher than the experimental level at predicted separation. 

6. Satisfactory estimates of the pressure recovery coefficient downstream of the 
predicted separation point were obtained by applying the constant effectiveness cri- 
terion. The best results were obtained when the predicted separation was the most 
premature. The agreement in this case was within about 12 percent of the experi- 
mental values of pressure recovery. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 13, 1973, 

501-24. 
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APPENDIX A 


HOT-WIRE EQUATIONS 

Two t 3 ^s of hot-wire probes were used to obtain the measurements in this study; 
namely, single wires normal to the flow and X- wires. The single wire probes were 
used to measure the mean velocity profiles and loi^tudinal component of turbulence in- 
tensity and the X-wires provided information from which the lateral component of tur- 
bulence intensity and Re 3 molds stress could be estimated. Since the derivation of the 
equations is presented in several publications, only the final form of the equations which 
were used in this study will be presented. 

Each wire was calibrated at the diffuser centerline to obtain the bridge voltage E 
as a function of the velocity U which were related according to the following general 
expression: 


E^ = (a + bu")(T^ - Tg) (Al) 

where a, b, n, and (T^ - Tg) are constants. The temperature difference (T^ - Tg) was 
nominally 110 K and the value of n in equation (Al) was 0. 5. Also during calibration a 
linearizer with zero velocity suppression was used to obtain an output signal which was 
directly proportional to U. The output signal from the linearizer was fed to an x-y 
recorder to obtain direct plots of the mean velocity profiles U/Ug(x) = f[y/R(x)] in the 
diffuser . 

The longitudinal component of turbulence intensity u’/Ug(x) was calculated from the 
linearized form of the hot-wire equation 



Ug(x) Ug(x) 



(A2) 


where represents the bridge voltage at zero velocity and Ve^ is the root- mean- 

square (rms) component of bridge voltage. All rms quantities were measured with a 
3-second time constant circuit. 

The signals from the X-wire probe were fed to a commercial signal analyzer and 
correlator which contained a sum and difference network. The output signals from the 
correlator were used to calculate the lateral component of turbulence intensity v'/Ug(x) 

and shear stress function uv/Ug(x). The expressions for these quantities as given in 
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reference 18 are 



The subscripts I and n in equations (A3) and (A4) refer to the two wires of the X- 
array and the angle <p {cp = 45® ) is the angle between the wires and the direction of the 
mean flow. A discussion of the effect of X-wire misalinement with the mean flow on the 
measxxred values of uv is given in appendix B. 
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APPENDIX B 


EFFECT OF X-WIRE MISALINEMENT ON MEASURED VALUES OF SHEAR 
STRESS AT THE DIFFUSER CENTERLINE 

In an axisymmetric flow field having a boundary layer which extends to the duct 
centerline, the shear stress term uv should be zero at the centerline and change sign 
in crossing the centerline. In the present study, the effect of X-wire misalinement with 
the mean flow direction in the vicinity of the diffuser centerline is believed to be the 
principal reason for observed nonzero values of uv. This alinement error can occur at 
stations 9 and 17 in the diffuser since the X-wire probe is inserted into the stream in a 
direction normal to the wall. The effect of the resulting misalinement at the centerline 
can be estimated for these two stations by accounting for the imequal angles of the two 
wires with the mean flow direction. In appendix A the hot-wire equations were based on 
equal angles between the wires and mean flow direction. If these angles are unequal, 
the effective velocities sensed by the wires are given simply as the vector sum of the 
velocity components as shown with the aid of the following sketches: 



Wire number I 


vf>'^ 

Wire number II 


Thus, 



{[(U + u) sin a + V cos a]' 


+ w 


,11/2 


(Bl) 


and 




+ u) sin ^ + V cos + w‘ 


1 


1/2 


(B2) 


where w is the fluctuating component of velocity perpendicular to u and v. In the 
foUowing derivation only wire I will be considered since the equations pertinent to 
wire II can be obtained by the same procedure. 

Equation (Bl) (dropping the subscript I) can be expanded in a binomial series to 

give 
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= (U sin Qf)*^ fl + — + — cot 0 !^ 
\ U U / 


(B3) 


where it is assumed that w « U. The equation for the hot wire relating the instanta- 
neous voltage E to the effective velocity V is 


E2=(aj*bjV”)(T^-T^) 


(B4) 


where (T^ - T^) is held constant. Substitution of equation (B3) into equation (B4) yields 


e2 = 


Lj + bj^(U sin ci!)*^^l + — + — cot O!^ ^(T. 


T ) 
w e' 


(B5) 


By analogy with equation (Al) 


E^ = [aj + bj(U sin - Tg) 


(B6) 


when U = 0, 


or 


e2 = e 2^, . aj(T„ - T,) 


- ®ref = ■ "^e^ 


(B7) 


If it is assumed that 


fl + — + — cot a\ «l + nf— + — cot o;) 
\ V U / \U U / 


(B8) 


substitution of equations (B8), (B7), and (B6) into equation (B5) yields 


2eE 


u V . 

= — + — cot a 




(B9) 
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Now consider wires I and II and define 


and 


ZejEi 






^®n®ii 




.n) 


U V . 

= — + — cot a 
U U 


= ^^--^cot;8 
U U 


Squaring and averaging equations (BIO) and (Bll) yields 


(BIO) 


(Bll) 


X2- 


4efE2 


" ®ref,I 


2 2 

= — + — cot“ a + 


cot a 


(B12) 


and 






n2/k2 _ p2 

” '®n ®ref , n 






U2 


2uv 

U2 


cot j3 


(B13) 


Equation (B13) can be subtracted from equation (B12) to give the following relation for 


w u2(x2 - Y^) - v^(cot^ a - cot^ /3) 
2U^(cot a + cot ^) 


(B14) 


2 

The value of v is determined by subtracting equation (Bll) from equation (BIO), 
squaring, and averaging. Then 


- 2XY 
(cot or + cot /3)^ 


(B15) 
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The quantities and in equations (B14) and (B15) can be expressed in terms of 
the measured voltages by means of equations (B12) and (B13). In order to express the 
product XY in terms of measured quantities it is necessary to resort to equations (BIO) 

and (Bll). From these equations the sum ^ej + and difference (^I - are 
formed and subtracted to give 


_ ^ (ej + ejj)^ - (ej - 



(B16) 


Substitution of equations (B15) and (B16) in equation (B14) leads to the following equa- 
tion for u^U^(x) in terms of the measured voltages and wire angles: 


uv 

U2(x) 


in^(cot a + cot /3) 


2^72 

®I^I 


2:pr2 

®n^n 




2 (cot a - cot )3) 

.2/ , . o\2 


e^E^ 


e^ 

®ii^n 




i®i®n 

(°i ■ (°i ■ °n)^. . 

‘{n- 



u 


Uc(x) 


(B17) 


Note that when oi = ^, equation (B17) reduces identically to equation (A4). 

When the X-wire probe is located at the centerline of the diffuser (stations 9 and 17) 

wire number I forms an angle of 45. 0° + 6. 5° =51.5° with the mean flow direction and 

wire number II forms an angle of 45. 0° - 6. 5° = 38. 5° with the mean flow direction. 

Using the data for these stations in conjimction with equation (B17) the calculated values 

of uv/U^(x) become +0. 0006 and -0. 002 for stations 9 and 17, respectively. The im- 
• ^ / 2 

corrected values of uwU^(x) for the above stations were -0. 0008 and -0. 006 as shown 
in figure 9. In both cases the angle correction tends to drive the value of uv towards 
zero as expected. 
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A corresponding expression for v'/U (x) can be obtained from equation (B15) and 

the expressions for X^, Y^, and XY given by equations (B12), (B13), and (B16), re- 
spectively. The equation for the lateral component of turbulence intensity becomes 


r 


V* 

Uc(x) 


®I^I 


®n^n 


ln^(cot a + cot j3)^ 





(®i ■ (®i ■ ®n)^_ 

2 (vj- 



0.5 


U 


U,(x) 


(B18) 


This equation is identical to equation (A3) except for the factor involving the angles. 
Dividing equation (B18) by equation (A3) gives the ratio of v’ for unequal wire angles 
to the value of v’ based on equal angles. This ratio is simply 2/[tan cp {cot a + cot j3)]. 
Substitution of ^ = 45° (appendix A) and the previous values for a and jS into this 
term indicates that at the centerline the corrected value of v' is only 3 percent less 
than the measured value. 
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